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ABSTRACT 

We explore the effects of small scale structure on the formation and equilibrium of dark matter 
halos in a universe dominated by vacuum energy. We present the results of a suite of four N-body 
simulations, two with a ACDM initial power spectrum and two with WDM-like spectra that suppress 
the early formation of small structures. All simulations are run into to far future when the universe 
is 64/i~^Gyr old, long enough for halos to essentially reach dynamical equilibrium. We quantify the 
importance of hierarchical merging on the halo mass accretion history, the substructure pop ulation, 
and th e equilibrium density profile. We modify the mass accretion history function of Wcchsler et al. I 
(|2002f ) by introducing a parameter, 7, that controls the rate of mass accretion, d\iiM/d\na cx a"''', 
and find that this form characterizes both hierarchical and monolithic formation. Subhalo decay rates 
are exponential in time with a much shorter time scale for WDM halos. At the end of the simulations, 
we find truncated Hernquist density profiles for halos in both the CDM and WDM cosmologies. There 
is a systematic shift to lower concentration for WDM halos, but both cosmologies lie on the same 
locus relating concentration and formation epoch. Because the form of the density profile remains 
unchanged, our results indicate that the equilibrium halo density profile is set independently of the 
halo formation process. 

Subject headings: cosmology:theory — large-scale structure of universe — dark matter 
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1. INTRODUCTION 

Our understanding of cosmology has consolidated sig- 
nificantly in recent years and we are now approaching 
percent-level estimate s of the most import ant cosmolog- 
ical parameters fe.g.. ISpergel et al~ll2006( l. Using these 
parameters, numerical simulations of large scale struc- 
ture provide excellent descriptions of the distribution and 
properties of dark matter halos. The next open question 
is to understand why non-linear structure takes the form 
that it does — a form that has been predicted by sim- 
ulations and confirmed, more or less, by observational 
data. In particular, we would like to know how and why 
dark matter halos attain a nearly universal form for their 
density profiles, as first described by Navarro, Frenk, & 
White 1997 (hereafter NEW). One aspect of this issue 
is understanding the importance of the method of mass 
accretion: How much does the final structure depend 
on accreting mass as virialized clumps as opposed to a 
continuum of diff use material and ho w effectively does 
violent relaxation (|Lynden-Bell I[l967l ) erase the memory 
of this accretion process. 

The basic process for the buildup of structure in 
our cold dark matter (CDM) dominated universe 
is the h ierarchica l merging of collapsed structures 
Csee. e.g..lPress fc S chechter 19 71 lAarseth et al. Ill979t 
iBlumenthal et al. 111984 : .Davis et aL I Il985f ). This pro- 
cess creates small halos early in the universe which 
merge with each other while accreting material from their 
surroundings, eventually creating the large cluster-sized 
structures of today through a "bottom up" process. 

^ Astronomy Department, University of Michigan, Ann Arbor, 
MI 48109 



For a time it was thought that light neutrinos 
might dominate the mass density, forcing galaxy for- 
mation to occu r thro ugh a "top down" process (e.g., 
iBond &: Szalav I Il983l ). In such a hot dark matter 
model, perturbations on small scales are washed out 
by free streaming, preventing the formation of the 
early low-mass seeds of hierarchical structure forma- 
tion. Dark matter halos still form, but the first ob- 
jects are large cluster-mas s halos (iZel'dovich I [l970l: 
iDoroshkevich fc Zel'dovichl Il975f ). While popular in 
the 1970's because simulations reproduced the outline 
of the cos mic web that surveys were .j ust beginning to 
map out (jThompson fc Gregorv I [l978l ). hot dark mat- 
ter models have been rul ed out based on observations o f 
the galaxy distribution (jWhite. Frenk. fc Davis 1 11983"). 
Tuning the free-streaming mass-scale leads to WDM 
possibilities which suppress density perturbations below 
some (typically dwarf galaxy sized) scale. The most 
immediate effect of this suppression is to reduce the 
number of small halos and subhalos existing in large 
halos. Several numerical stud ies of WDM cosmolo- 
gies have been carried out (e.g.. lEvrard fc Crone | [l992 



Avila- Reese et al . 1120011: 'Bode. Ostriker. fc Turok 



2001 



Coli'n Avila-Ree se. fcValcnzucla I I2000I: IWhite fc Croft I 
20001 iKnebe et al. 112002 . [2003), mostly in an attempt 



to explain the apparent lack of substructure in our local 
group as compared to predictions from ACDM simula- 
tions. 

Although CDM and hierarchical merging has emerged 
as the standard paradigm for structure formation, and 
is the key ingredient in setting the distribution of dark 
matter halos, it is still uncertain how important this pro- 
cess of mass accretion is in setting the internal proper- 
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ties of dark matter halos. While some studies indicate 
that the accretion of substructures plays a significant 
role in sett ing the inner slope of the r adial halo den- 
sity profile ()Ma &: Bovlan-Kolchin |[2004[ ). simulations of 
monolithic collapse events in WDM scenarios produce 
halos with global properties unch anged relative to similar 
structures in CDM cosmologies (|Evrard fc Crone Ill992l : 
iMoore et al. 1 11999| ). By truncating the initial power 
spectrum in an otherwise CDM simulation at some scale 
kc (with corresponding mass scale Mc = A/Sn'^k'^p), one 
can mimic WDM cosmologies and test the importance 
of hierarchical growth for establishing halo structure. 
Compared to the CDM model, where mass is continu- 
ally accreted in dense clumps, WDM cosmologies have 
an early period of monolithic collapse, where large halos 
form out of a smooth background and relax with many 
fewer disruptions due to merger events. If the halo is in 
a dense enough region, this initial collapse is followed by 
hiera rchical accretion. Previous simulations of this pro- 
cess (Moore et al. "1999 HBode. Ostriker. fc Turok II2001I : 
lAv ila-Recse et al. 2001]) have shown that the resulting 
density profile is virtually unchanged for halos well above 
the truncation scale. 

Halos in a WDM cosmology are effectively a re-scaled 
versions of the first halos expected in a CDM cosmol- 
ogy. Most physical CDM candidates (i.e., SUSY-LSP's) 
have some intrinsic velocity dispersion, washing out per- 
turbations on very small scales (much smaller than we 
are able to simulate on cosmological scales), effectively 
trun cating the power s pectrum at som e very large k 
(see iDiemand. Moore, fc Stadcl 2005; Gao et al. 1 12005| 
for discussions of such simulations). In this manner, 
studying WDM cosmologies can provide clues to the ear- 
liest collapse of CDM structures. 

The discovery of the accelerated expansion of the 
universe (Perlmutter et al. 1998), potentially caused 
by a non-zero cosmological constant (the ACDM 
model), provides us with a mechanism for studying 
the true equilibrium structure of cosmic halos. Previ- 
ous w ork f Adams fc Laughlin "1 9971 : iNagamine fc Loeb I 
120031 Busha ct al. 2003, 2005) has shown that, in the 
relatively near cosmic future of a ACDM universe, struc- 
ture formation and halo growth will come to a rapid end. 
This happens around a ~ 3, when Q\ « 1. At this time, 
exponential deSitter expansion causes mergers and accre- 
tion to stop and pushes existing halos further and further 
away from each other. From a ~ 3 beyond, no matter 
is left for halos to accrete, and they relax toward their 
asymptotic equilibrium state in effective isolation. These 
late-time halos experience no disruption from mergers or 
other events that frustrate the equilibrium, in contrast 
to the situation at the present epoch. 

To gain insight into the question of the origin of the in- 
ternal structure of dark matter halos, we have performed 
a set of simulations of cosmic structure, using both a 
standard ACDM and truncated WDM-like power spec- 
tra. Simulations are run into the far future, allowing 
halos to relax toward equilibrium configurations. Our 
numerical simulations are described in §2. In §3, we com- 
pare the differences in the halo distribution for the two 
cosmologies, concentrating on the mass function, evolu- 
tion of the power spectrum, and the formation of WDM 
halos with mass well below the truncation scale. In §4 
we compare the internal structure of the halos in the 



TABLE 1 
Simulation parameters: 



L[h- '-Mpc] 


Mplh~'MQ] 


ep[/i~-^Mpc] Np 


200 


3.97 X IQi" 


0.04 256^ 


50 


6.20 X 10* 


0.01 256^ 



two cosmologies, including the mass accretion histories 
(MAHs), the distribution of substructure, and the halo 
density profiles. Results and their implications are sum- 
marized and discussed in §5. 

2. SIMULATIONS 

We simulate the formation of dark matter halos in 
A-dominated CDM and WDM-like cosmologies with a 
suite of dark matter N-body simulations using the pub- 
hcly available TreePM code Gadget 2.0 (Springel 2005). 
We use two ACDM and WDM simulations with different 
mass resolutions for a total of four large-scale cosmolog- 
ical simulations. The simulation pairs at each mass res- 
olution were created with the same initial phases so that 
the the large-scale environment would be unchanged. All 
simulations model a patch of space in a flat universe 
with current matter density f2m = 0.3, vacuum den- 
sity fljy — 0.7, Hubble parameter Hq — 70km s^^Mpc^^ 
and power spectrum normalization erg = 0.9, values con- 
sistent with the first year release of WMAP me asure- 
ments of the CMB power spectrum (jSpergel et al. | [2003') . 
The lower resolution simulations model periodic cubes 
of side length L = 200/i^iMpc containing Np = 256^ 
particles of mass 3.97 x IO^'^H^^Mq and gravitational 
softening length Cp = 40/i~^kpc. The higher resolution 
simulations use cubes of side length L — 50/i~^Mpc, 
containing 256^ particles with mass 6.20 x 10*/i~^Mq 
and softening length ep = lO/i^^kpc. The softening 
scales quoted here correspond to their values at the 
present epoch and are held constant in comoving space 
for a < 1, but become physical lengths for a > 1 to 
prevent structures from being over-softened due to the 
exponential increase in a during the deSitter phase. Ta- 
ble [1] lists these simulation parameters. All simulations 
were started at redshift z = 19 (scale factor a = 0.05) 
and were run into the far future, a = 100. Although 
the starting redshift is somewhat late, it is consistent 
with the analysis by iPower et al. I ([2003) for a simula- 
tion of our resolution, and should present no problem for 
the late-time results we are primarily concerned with. 
We store a total of 300 outputs equally spaced in log(a) 
for each simulation. At a = 100 the universe is about 
64/i~^Gyr old and structure formati on has preceded to 
compl etion in a A CDM cosmology (Bus ha et alTI l2003t 
.Nagamine fc Loebl[200l ). The simulations were run on 
16 nodes of a dual-Opteron Beowulf cluster at the Uni- 
versity of Michigan. 

The initial power s pectrum for our CDM simu lations 
was set using cmbfast (S eljak fc Z aldarriaga' 1996) in ac- 
cordance with WMAP year 1 data (S pcr g el et a l. 2003). 
Our second, WDM-like, simulation used a truncated ver- 
sion of this initial power spectrum, 

Pt(fc) = CPo(fc)e-('=/^^)', (1) 

where kc is the truncation scale and C a normalization 
coefficient that allows us to set ag. Figure [1] plots our 
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Fig. 1. — The input power spectra for the ACDM (sold hne) 
and WDM (dashed line) simulations. The ACDM spectrum is cal- 
culated using CMBFAST . The WDM model adds an exponential 
cutoff to the ACDM model at a mass scale of 8.09 X lO^^h-'^MQ. 
Both spectra were normalized so that erg = 0.9. The light gray 
boxes represent the range between the fundamental and Nyquist 
frequencies for the small (top) and large (bottom) volume runs. 
The hatched-line boxes show the range between the Nyquist and 
softening frequencies for a < 1. For a > 1, the softening is constant 
in the physical frame and its corresponding wavenumber grows with 
fc DC a in this comoving representation. 

input spectra, with Pq as the sohd line and Pt as the 
dashed hne. For our simulations, we choose a truncation 
scale of kc = 0.511/iMpc-"'^, which corresponds to a mass 
scale of Mc = 8.09 x IO^/i-IMq (2037 particles for our 
lower resolution simulations and 130,396 at the higher 
resolution). This truncated spectrum was re- normalized 
to CTg = 0.9 so that high mass halo abundances would 
be similar. As noted earlier, phase information was re- 
tained for each resolution pair, resulting in similar large 
scale structures (Figure [2]) and allowing us to identify 
corresponding halos in the two cosmologies. 

Equation ([T]) represents a transfer function 
that differs from the standard transfer function 
for a WDM cosmology. Generally, the mass of 
the WDM particle, mw sets a free streaming 
length, Rf = 0.2(17w / i^)i/3 (mw/lkeV)-4/3Mpc 
(jSommer-Larsen fc Dolgovl l2001h . which approxi- 
mates the WDM pow er spectrum through the relation 
(|Bardeen et al. Ill986f ) 

PwDM{k) = exp [-kRf - [kRff] PcDM{k). (2) 

This spectrum has a slightly more gradual cutoff than 
equation It should be noted that we refer to our 
truncated models as "WDM" cosmologies even through 
they were not created using this transfer function. For 
reference, our truncated model most closely approxi- 
mates a WDM cosmology with mwDM = O.lSkeV and 
Rf = 1.6Mpc. 

The evolution of the resulting density fields of the sim- 
ulations are shown in Figures [5] and [31 Figure [5] shows 
a slice the comoving density field from our smaller vol- 
ume CDM and WDM simulations at a = 0.3,1,3, and 
100. The differences between these two models are strik- 
ing, especially at a = 0.3, where the CDM cosmology 
exhibits a well-formed web with an abundance of small 
halos. The WDM cosmology, in contrast, has a mostly 
uniform density, with only one visible halo and a hand- 
ful of weak filamentary structures making up the cosmic 



web. A clear cosmic web does rapidly develop in the 
WDM simulation, however, and by a = 1 similar large 
scale structures are present in both cosmologies, even 
though there is a strong suppression of small halos in 
the WDM filaments. By a = 3, the large scale density 
field is set and undergoes little evolution from a = 3 to 
100. Once the cosmological constant becomes dominant 
the growth function saturates, ending halo formation and 
freezing the comoving web. Halos continue to contract 
in this comoving picture, and by a = 100 they consist of 
small, tightly bound knots along and at the intersection 
of filaments. 

Figure [3] shows the evolution of the density field in 
a fixed physical region. Here, we focus on the evo- 
lution of a particular CDM halo with mass M200 = 
5.38 X lO^^h-^MQ and its counterpart in the WDM cos- 
mology (see §3.3). While much of the late time growth 
is identical (such as the major merger around a = 1), 
the initial formation processes differ substantially. At 
a = 0.3, there are many low mass CDM progenitor ha- 
los present, while the WDM halo looks like a weak (un- 
collapsed) perturbation in an otherwise smooth back- 
ground. This difference is manifest throughout all plot- 
ted epochs by the persistent lack of substructure in the 
WDM halo at the three later epochs. 

One way to quantify the expected suppression of hi- 
erarchical buildup is to look at the critical mass scale, 
(a) , where one expects a perturbation to go non-linear 
and collapse, defined through the relation 

<^[MM\ = (3) 

Here, a{M) = {2tt)-^ P(k)W'^{M,k)d^k, W is the 
Fourier transformation of the top-hat window function, 
D{a) is the linear growth function, and 6c = 1.686 is the 
linearly extrapolated crite rion for collapse of an over- 
dense perturbation (Press fc Schechter I Il97l : [Peebles I 
[19801) . The factorD(a) can be calculated numerically 
from the expression 

and is normalized such that D{1) — 1. The left panel 
of Figure [H shows the amplitude of 2ct(A/) perturbations 
as a function of M at the present epoch. The horizon- 
tal dotted line shows the critical scale, 5c = 1.686. The 
right panel of Figure [H shows the evolution of M2cr (a) 
— the mass of a 2a perturbation that goes nonlinear as 
a function of cosmic time. In both cases, the solid line 
represents the full power spectrum model and the dashed 
line represents the truncated power spectrum model. We 
consider 2a perturbations because, in our WDM cosmol- 
ogy, la perturbations do not collapse until a — 1.33 and 
are prevented from growing substantially by the cosmo- 
logical constant. The normalization erg — 0.9 boosts the 
power at large mass (M > IO^/i-^Mq) in the WDM 
model relative to the CDM case and causes the largest 
perturbations to collapse slightly earlier and with larger 
asymptotic masses in the WDM model. We also expect 
no structures in the WDM cosmology for a < 0.4, agree- 
ing with Figure [31 in which the most massive halo of the 
simulation has yet to form a coherent structure prior to 
this epoch. 
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Fig. 2. — The density fields of comoving slices of the cosmic web at o = 0.3, 1,3, and 100 (columns, left to right) of a ACDM (top) and 
WDM (bottom) cosmologies from our small volume simulations. The side length for each image is 35/i~^Mpc, with thickness 5/i^^Mpc. 
The grey-scale is proportional to log(p/p). 
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Fig. 3. — The density field around a large halo at a = 0.3,1,3, and 100 (columns, left to right) from our small volume simulations. 
The top row shows the largest halo from the ACDM cosmology (M200 = 5.38 X 10^*/i~^Mq at a = 100) and the bottom row shows the 
corresponding WDM halo (M200 = 5-24 X 10^*/i~^Mq at a = 100). The side length for each image is 12h~^Mpc, with thickness 6/i~^Mpc 
in physical units. The grey-scale is proportional to log(p/pc). 



Dark matter halos in our simulations are identified us- 
ing a standard friends-of-friends (FOF) algorithm with 
linking length 0.15 times the inter-particle spacing. Halo 
centers are identified as the most bound particle of the 
resulting group. As a mass measure capable of spanning 
both early and late times we use M200, the mass of all 
particles inside a sphere of radius r2oo with over-density 
200 times the critical density, pc- In previous work, we 
found the M20Q provides a good proxy for the asymptotic 
halo mass and is roughly half the value of the ul timate 
halo mass, M200 ^ Q.^Muaio (|Busha et al. 1 120051) . The 
halo velocity is defined to be the center of mass velocity 
of all particles within r2oo • 

Subhalos were ide ntified using the subfind routine 
(|Springel et al. Il2001h . This routine works on top of an 



FOF group and identifies density maxima within halos 
from an SPH smoothing kernel that uses the distance 
to the 32"''-nearest neighbor to obtain local density esti- 
mates. Subhalos are then selected as locally overdense re- 
gions containing at least 20 bound particles. The largest 
subhalo identified by subfind is actually the host halo 
of the FOF group. At a = 100, this host halo should 
correspond to the actual equilibrated halo as defined in 
iBusha et al.1 ([2005,), minus any locally bound subhalos. 
A comparison between these mass estimates shows good 
agreement, with the masses agreeing to within 1%. 

Throughout this paper, we use several scale radii. Halo 
sizes are defined using r2oo (see above) and Vhaio, the 
spherical radius containing all bound particles (which is 
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Fig. 4. — Left: Rare (2(t) perturbation amplitude as a function 
models. The vertical lines show the mass scales for 100 particles at 
of scale factor for the ACDM (solid line) and WDM (dashed line) 
particles at the two resolutions. 



of mass at a = 1 for the ACDM (solid line) and WDM (dashed line) 
the two resolutions. Right: Characteristic collapsed mass as a function 
cosmologies. The solid horizontal lines show the mass scales for 100 



only defined for a ^ 5, see lBusha et al. Il2005f) . Addition- 
ally, we fit NFW and Hernquist density profiles to our 
halos, 



Pnfw - —, — 77- — -, — r^, (5) 

r/rs{l + r/rsY 

which adds the scales (the radius where the best fit 
NFW profile has logarithmic slope —2) and Tc (where 
the best fit Hernquist profile has slope —2.5). Generally, 
we find that rhaio — 4.6r2oo and rs = 0.4rc (see §4.3, 
iBusha et al.1 12005). 

3. COMPARISONS OF THE DARK MATTER 
DISTRIBUTION 

In this section we compare properties of the distribu- 
tion of dark matter halos in our simulations, including 
the evolution of the power spectrum, the halo mass func- 
tion, the correspondence between CDM and WDM halos, 
and the formation of WDM halos below the truncation 
scale. 

3.1. Evolution of the Power Spectrum 

Figure [5] shows the evolution of the power spectrum for 
all four of our simulations at the epochs a = 0.3, 1, 3, and 
100. The solid lines represent the CDM cosmologies and 
the dashed lines the WDM spectrum. The power spec- 
trum is shown in dimensionless units, A^(fc) oc k^P{k). 
Power spectra for the large and small volume simulations 
are combined, allowing us to probe a larger range in k. 
The arrows represent the softening lengths for the large 
and small volume simulations at the plotted epoch. The 
collapse of non-linear structure creates substantial power 
beyond the Nyquist frequency, which we measure using 
the tiling method of I Jenkins et al.l ()1998[ ). The spectra 
are plotted from the fundamental simulation frequency 
out to a wavcnumber where the shot noise of a Poisson 
distribution of particles becomes comparable to the mea- 
sured power. 



At a = 0.3, when non-linear structure formation is in 
its early stages, the power spectrum of the WDM model 
is heavily truncated above fee. By the present epoch, 
much of this suppression has disappeared due to power 
transfer from collapsing structures, and the WDM cos- 
mology matches the CDM model almost perfectly at low 
and intermediate wavenumber, up to an order of mag- 
nitude above kc- Relatively little happens to the power 
spectrum beyond a = 1. As noted earlier, the dominance 
of A halts the growth of structure beyond a ^ i and 
causes the power spectrum to freeze after only a modest 
amount of additional evolution. Since Figure [5] plots the 
power spectrum in comoving space, beyond a = 1 the 
expansion of the universe transfers power to larger scales 
with no real change in the shape of the spectrum. For 
most of the measurable range, 0.1 < k < lOOhMpc"^, 
the power spectrum is well characterized by the expected 
power law, cx k^, caused by the freezing of struc- 
ture in an expanding universe. Collapse, however, has 
not managed to elevate the WDM power spectrum up 
to that of the CDM model at all scales by a = 3, and a 
small suppression of power at the highest wavenumbers 
persists at a = 100. This suppression is due to the lack 
of low mass halos. If we measure the halo power the 
spectrum using only halos with M200 > Mc we recover 
spectra that are identical at the few percent level, which 
indicates that the distribution of halos on scales greater 
than Mc is statistically unchanged. 

The evolution of the power spectru m in WDM cos- 
molog ies was studied p reviously by IWhite fc Croft I 
(|2000l) and iKnebe et al. I ()2002[ ). These studies present 
results similar to ours. By a ~ 0.5, non-linear collapse 
has boosted the truncated portion of the power spec- 
trum substantially, agreeing with a ACDM model for 
k < lO/iMpc""'^ with only a slight suppression for larger 
k. 



3.2. Mass Function 
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Figure [S] shows the mass function for our large volume 
CDM and WDM simulations at a = 1 and 100. Errors 
are calculated assuming Poisson statistics and the verti- 
cal dotted line is the truncation scale, Mc- Here, the dark 
lines represent the ACDM model and the light lines the 
WDM version. For comparison , the figure also shows the 
Jenkins Mass Functions fJMF. lJenkins et al. I ()2001f )) as 
dashed lines. The JMF is defined via 

dnjMF{M, a) p dlno'^"'^(M, a) 

■ — ^ ~ X 



Fig. 5. — The evolution of the power spectrum for halos in CDM (solid lines) and WDM (dashed lines) cosmologies. The spectra are 
shown at a = 0.3, 1, 3, and 100. Power spectra from the small and large volume realizations are laid on top of each other. The dotted 
vertical lines represent the truncation scale, kc, and the arrows show the softening scales for the large and small volume runs. The spectra 
are plotted from the fundamental frequency of the simulation volume to the frequency where the shot noise becomes comparable to the 
measured power, generally ~ O-lkg^ft- 

one particular realization. The use for FOF masses in 
measuring the mass function does create difficulties in 
the far future because the FOF groupfinder identifies 
isosurfaces relative to the background matter density, 
which is dropping rapidly from the exponential expan- 
sion. Compared with the critical density, the physi- 
cal density isosurface identified by a FOF groupfinder 
is pfof = [ia^^Tn,o)/i^Tn,o + a^^A,o)]b~^Pcrit- By a = 3, 
Pfof — 3.5pcrit, which includes a substantial amount 
of material outside the virialized region of a halo that 
is unbound and bei ng pulled away by the Hubble flow 
(|Busha et al. Il200"5[ ). At late epochs, Mpop is a poor 
proxy for the actual (virialized) mass of a halo. The 
M200 mass function, however, doe s not suffer from this 
defect. The a = 1 result agrees with' Evrard et ahl ()2002f ) 
to ~ 20% in number, approximately the quoted statisti- 
cal accuracy. The parameters are slightly off, however, 
due to insufficient mass resolution of the simulation. The 
a = 100 mass function does provide a substantially bet- 
ter agreement than the FOF mass function, fitting to 
within about 2a at all masses. 

The mass function for the WDM cosmology exhibits a 
striking suppression for low masses, beginning slightly 
above the the truncation scale, Mc (dotted vertical 
line). An unexpected upturn appears in the WDM 
mod el at the lowest masses ( ^ 50 par ticles ) for FOF ha- 
los. iBode. Ostriker. fc Turok I (|2001h and iKnebe et al~l 
(I2OOI " claim that this behavior is a result of physical 



d\n{M) "M dlnM 

exp[-\lna-\M,a)+B\'], (7) 
where A, B, and e are fitting parameters. Two mass val- 
ues are used in this plot to compare with published 
JMF parameters: FOF ma sses with a lin king length b — 
0.164 (upper curves, Jenkin s^t al. Il200l[ ) and M200 from 
a spherical ov erdensity groupfinder (SO, lower curves, 
lEvrard et al. |[2002i) . The fitting parameters are listed in 
Table H For the a = 100 spherical ove rdensity JMF 
we us ed the O^/ = parameters from lEvrard et al. I 
(|2002f) . Ahhough not shown. The She th fc Tormen mass 
funct ion with published parameters (jSheth fc Tormeiil 
|1999| ) agrees quite well with the JMF FOF mass func- 
tion at all epochs. While the agreement between FOF 
masses in the CDM cosmology and equation ([7]) is good 
for a = 1, a substantial mass excess is present in 
the mass range 10" < Mfof < 5 x IO^'^/i^^Mq at 
a = 100. Additional simulations of this cosmology con- 
firm that this excess is significant and not a result of 
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Fig. 6. — The mass function for the CDM and WDM cosmologies (dark and hght curves) at a = 1 and 100 for our larger volume 
simulations. The upper curves show the FOF{0.164) mass function and the lower curves the SO(200) mass function. The FOF(0.164) 
function has been offset vertically for clarity. The vertical dotted line shows the truncation scale. Error bars assume Poisson statistics. 



TABLE 2 
Mass Function Parameters: 



Mass Function 


A 


B 


e 


JMF - FOF(0.164) 


0.301 


0.64 


3.88 


JMF - SO(200, a = 1) 


0.220 


0.73 


3.86 


JMF - SO{200, a = 100) 


0.199 


0.76 


3.90 



halos forming through fragmentation from Jeans insta- 
bihty. However, such an upturn is not present in the 
SO mass function and we show evidence in the ap- 
pendix that it is actually a numerical artifact of the 
FOF groupfinder. We have also calculated JMF and 
ST fits for our WDM cosmology, which do not fit as 
well as in the ACDM cosmology. Both the JMF and ST 
mass functions strongly over-predict the abundance of 
halos with masses M < Mc- The poor fit in the range 
IO^^/j-IMq < M < 8 X lO^/i-iM© should not be sur- 
prising because equation ^ was motivated by a pertur- 
bation collapse threshold, similar to the Press-Schechter 
(1974) model, which uses spherical collapse to determine 
a collapse epoch. As we discuss in the next section, 
WDM halos with mass below Mc form out of larger mass 
perturbations that do not follow this model, at odds with 
the assumptions of Press-Schechter. 

3.3. Halo Correspondence 

Because we used the same phases in constructing the 
initial conditions, we can cross-match halos in the CDM 
and WDM cosmologies using a Lagrangian scheme. We 
select a FOF halo in one cosmology (usually ACDM ) 
and identify for the largest halo in the other cosmology 
containing at least 50% of particles of the selected halo. 
When starting with a ACDM halo, we don't allow any 
two WDM halos to be identified with the same ACDM 
halo. This simple method is robust for massive systems, 
and corresponding halos are found for 98% of all halos 
with M200 > Mc = 8.08 X 10"/i-1Mq. The "missing" 
halos are lost because the smoothing of the power spec- 



trum causes distinct halos in the CDM cosmology to form 
as single halos in the WDM run. Our requirement that 
each WDM halo have only a single CDM counterpart pre- 
vents all but the most massive of these CDM halos from 
having a WDM counterpart. If we relax this correspon- 
dence requirement, all CDM halos with MpoF > Mc 
have WDM counterparts. Figure [7| shows the masses 
of corresponding halos in the WDM and ACDM mod- 
els. At higher masses (M200 > 2Mc) the masses are 
similar, M2oo,cdm ~ M2oo^wdm with about a 6% scat- 
ter. As the mass falls below the truncation scale, Mc, 
the WDM halos become less and less massive relative 
to their CDM counterparts, eventually disappearing al- 
together. There are a few extreme outliers from the 
general relationship. Halos with a low Mwdm /Mcdm 
(the three halos with M2oo,cdm ~ 10^'^h~^MQ and 
M2oo,WDM 2 X 10i:^/i"1Mq) are CDM halos in fila- 
ments that never completely collapse in the WDM cos- 
mology. Most of their particles are located in a spray 
along the filament, but there is a small WDM halo in 
the filament with ^ 90% of its members in the CDM 
halo. The halo with the M2oq,wdm ^ M2w,cdm is a 
rare occurrence where many small CDM halos 10) 
had merged into a single, much more massive halo in the 
WDM cosmology. 

Also plotted in Figure [7] is the completeness function 
for identifying corresponding halos. Virtually all ACDM 
halos with M200 > Mc have corresponding WDM ha- 
los, but the completeness function drops very rapidly 
for M200 < Mc, the mass range where halo formation 
is strongly suppressed. 

3.4. Sub- Truncation Scale Halos In WDM 

One surprising observation from the WDM simulations 
is that many halos form with mass scale M200 well be- 
low the truncation scale, Mc- Although physical halos 
with M200 < Mc = 8.08 X W^H'^Mq account for only a 
few percent of all collapsed mass, the se halos repr e sent a 
substantial population by number. Knebe et al. I ()2003f ) 
have shown that these halos appear in well-defined fila- 
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Fig. 7. — A comparison of the Af200 values for corresponding 
halos in the ACDM and WDM runs at a = 100. The dashed 
vertical line denotes the truncation scale. Halos with M200 3> Mc 
have roughly the same mass with a scatter of about 6%. The light 
line shows the probability for finding a WDM halo corresponding 
to a ACDM halo of a given mass. 

ments an d argue that they are th e result of Jeans insta- 
bility (Bi nney fc Tremaine Ill98"7l ). Our results confirm 
that these halos form in filamentary structures, and we 
find that many of these halos have CDM counterparts, 
indicating that these low mass WDM halos are not nu- 
merical artifacts. 

In order to understand how these sub-truncation scale 
halos form, we identify particles in small halos (A/200 < 
Mc) in our small volume simulations at a = 100 and trace 
them back to the initial conditions. We then locate the 
density peak nearest the center of mass of this distribu- 
tion and measure the mass of the perturbation by find- 
ing the radius where spherical overdensity drops below 
ScD{a). For all halos in the WDM simulation, the per- 
turbation mass is in the range (0.2 — 2) x lO^'^/i^^Mo = 
0.25 — 2.5Mc. In particular, no halos form out of pertur- 
bations substantially smaller than our truncation scale, 
but sub-truncation scale halos do form out of large per- 
turbations that do not fully collapse. The CDM sim- 
ulation, by comparison, has small halos forming out of 
perturbations with masses anywhere between the final 
halo mass, Mhaio, and 2 x 10^'*/i~^Mq, with the bulk of 
the halos (80%) forming from perturbations with mass 
less than 4 times Af/ja;o(100). The distribution of per- 
turbations with initial masses greater than Mc, how- 
ever, is almost identical to that for the WDM cosmology. 
Low mass WDM halos are simply large perturbations 
that do not collapse completely, and behave the same 
as in a CDM cosmology. S uch a trend was noted by 
iKatz. Quinn. fc Gelb I |l993f ). and our WDM halos seem 
to be an extreme case of the tendency for halos to form 
objects with a substantially different mass than their ini- 
tial spherical collapse prediction. 

4. COMPARISON OF HALO PROPERTIES 

This section compares properties of individual halos in 
our CDM and WDM cosmologies. In particular, we focus 
on the mass accretion histories, subhalo abundance, and 
density profiles. In spite of some substantial and funda- 
mental differences in the first two properties, the form 



of the radial density profiles is unchanged between the 
CDM and WDM cosmologies. Furthermore, halo concen- 
trations follow the same relation with formation epoch in 
both models. 

4.1. Mass Accretion Histories 

One property for which we expect a clear difference 
between CDM and WDM halos is the halo mass accre- 
tion history (MAH). For the WDM model, the reduced 
merger activity should result in a smoother MAH since 
mass is primarily accreted in the form of diffuse material. 
The suppression of power at large k also alters the charac- 
teristic collapse mass at low mass/early times (Figure IH 
right panel). The implication is that halos will form later 
and more rapidly (in the sense of a larger din AZ/dlna) 
in a WDM cosmology. 

The MAHs are measured using a halo's most massive 
progenitor, where a progenitor is any halo at a preceding 
output in which at least 50% of the FOF particles end up 
in the subsequent halo. When comparing halos between 
the two runs, we first select CDM halos from a given mass 
range and then select either their WDM counterparts (as 
in Figure [7l ignoring CDM halos that have no match) or 
WDM halos from the same mass ranges. These selection 
methods are nearly degenerate for Af2oo > Afc, as shown 
in Figure [71 At lower CDM halo masses A/2oo,wdm is 
strongly suppressed, so the selection methods differ sub- 
stantially. Generally, we prefer to consider halos of sim- 
ilar final masses so that we do not have to worry about 
halos with no counterpart. 

Figure [8] shows MAHs for three halos in the full run 
(solid curves) and their matched halos in the truncated 
run (dashed curves). These individual halos were se- 
lected from our smaller volume run to have A/200 S> 
Mc (dark curves), A'f2oo ^ Mc (medium curves), and 
M200 < Mc (light curves) at a = 100 in the CDM cosmol- 
ogy. As expected, the halos in the truncated model have 
smoother MAHs during the initial halo growth phase 
and form slightly later, with these effects becoming more 
pronounced for smaller halos. The overall shape of the 
MAHs for the halos with A/200 ^ Mc is remarkably 
similar. While the most massive halo in the truncated 
model has no progenitor before a ~ 0.2, it grows quickly 
and catches up with the CDM halo by a = 0.7. Af- 
terwards, the two halos evolve almost identically, even 
undergoing the same major mergers around a = 1 and 2. 
These mergers happen slightly earlier in the WDM cos- 
mology due to the increase in power at this scale from 
our normalization = 0.9, which effectively starts the 
"cosmic clock" for these large halos at a later time. In 
contrast with the these late-time similarities, the MAHs 
at a < 0.7 are substantially different. In the WDM 
scenario, a large mass perturbation will collapse more 
or less as a unit, as soon as it goes non-linear. This 
collapse creates a phase of smooth, rapid mass growth 
with din A/wDM/'^liia ^ din McBM/dlaa. Once the 
halo reaches A/200, wdm ~ A/c, the halo is the approxi- 
mate size of its counterpart in the CDM cosmology and 
d In A/200/d In a drops to match the rate of the CDM 
halo. The WDM halo then begins to accrete mass as 
already-collapsed clumps in a quasi-hierarchical fashion. 
The CDM halos with masses A/200 < A/c, in contrast, 
have corresponding WDM halos that never accrete mass 
in virialized clumps. The MAHs of such WDM halos are 
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Fig. 8. — The accretion history for three individual halos in 
the full run (solid curves) and their corresponding halos in the 
truncated run (dotted curves). The halos were selected with 
A^200 3> Mc (dark curves), M200 ~ Afc (medium curves), and 
Af200 A^c (light curves) in the run with the full power spec- 
trum. The dotted horizontal line represents Mc. 

much smoother and accrete the bulk of their mass in a 
single period of rapid accretion. 

Figure [9] shows ensemble average MAHs for halos from 
our large volume simulations. Here, the diamonds rep- 
resent CDM halos and the crosses are WDM halos. Ha- 
los from different mass ranges have been offset in time 
to make the figure more readable. Both the CDM and 
WDM halos are selected to lie in the mass ranges M200 = 
(2 - 4) X 10i3/i-iMo,(0.5 - 1.3) x IQi^/i-iM©, and > 
4 X IQi^'/i^^Mo at a = 100. Figure [9] shows many of the 
same trends observed in Figure [51 For the most mas- 
sive halos, the two cosmologies again show mass equal- 
ity around M2oo,cdm « -M200,wdm ~ Mc, with the 
WDM halos accreting mass significantly faster before this 
time. For M200 < Mc, WDM halos form later and more 
rapidly than CDM halos with similar mass. Also shown 
in the bottom panel of this figure are the growth rates, 
d\n{M2ao) /d\n{a), of the halos and their fits to equation 
®. 

Wechsler et al. (2002, hereafter W02) proposed a fit- 
ting formula for the MAH of a halo up to the present 
epoch of the form 



(8) 



The free parameter S in this equation is used only in 
defining ac^wo2, the creation epoch for the halo, when 
d\nM/d\na = S. We choose to follow their conven- 
tion and adopt S — 2. This formula is fit to our CDM 
MAHs over the range a = 0.2 — 1.0 (but continued out 
to a = 100) and plotted as the dotted curves in Figure [S) 
In general, the fit is good for both CDM and WDM (not 
shown) models for a < 1, but overestimates halo masses 
by a factor of 2 in the CDM run and more than an order 
of magnitude for the WDM run at late times. If the fit is 
calculated for the full range, 0.2 < a < 100, the late time 
asymptote is correct, but d In M2oo/'^ In a is substantially 
lower than observed for either cosmology at all epochs. 
This behavior is probably an indirect result of the co- 
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Fig. 9. — Top Panel: Average MAHs for CDM (diamonds) 
and WDM halos (plus symbols) from our large volume simula- 
tions. The halos are selected from the mass ranges M200 = 
(1 - 4) X 10", (0.5 - 1.3) X IQl*, and > 4 X 10"/i-lMpc. The 
curves are offset in a by a factor of 2 (intermediate mass range) 
and 4 (low mass range) to make them easier to distinguish, but 
the short vertical lines represent a = 1 for all mass ranges. The 
dotted curves are fits to equations JSjl for the CDM halos and 
the solid and dashed curves are fits to equation (|9j for the CDM 
and WDM cosmologies. Bottom Panel: The mass growth rates, 
dln(M2oo)/'i\n{a), of the halos and plotted above and their fits to 
equation (|9)l. Curves are offset in both the horizontal and vertical 
directions to make them more distinguishable. 

incidence problem — the surprising observation that we 
live during the relatively short epoch where ~ ^A- 
The fit works well for a < 1, even when calculated for 
just a fraction of the region and then extrapolated. Al- 
though equation ([S]) was created for halos in a ACDM 
universe, a = 1 is not much later than the equality epoch, 
Qeq — 0.75, when f2,n — r^A- Once A becomes the dom- 
inant component of the universe, the gr owth function 
quickly saturates and halos cease to grow (|Busha et al. I 
2005) . Consequently, equation ^ approaches its asymp- 
tote much more slowly than halos feeling the full effects 
of a dominant cosmological constant. To capture the full 
histories in both cosmologies, we propose a generalization 
of the form 

^/Qe-('>c/ao)^f ((ao/a)-'-!)^ 



M{a) 



(9) 



where 7 is the rate index which sets the mass growth 
rate through the relation d\nM/d\na — S{ac/a)^. This 
variable is introduced such that Oc is still the epoch where 
M{a) has a logarithmic slope of S, but 7 sets how quickly 
a halo grows and asymptotes to its equilibrium mass. 
When 7 = 1 this generalized form reduces to equation 
([8]). We expect to recover 7 > 1, which corresponds to 
more rapid formation and faster asymptote behavior (see 
Figure El) 

The fits to equation ^ in Figure [9] (solid and dashed 
curves for CDM and WDM halos) closely follow the mea- 
sured MAHs. They provide good agreement over all 
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Fig. 10. — Comparison of the formation epochs defined by W02 
(equation [8]) and our generalization (equation [9]). The left panel 
represent halos from the CDM run and the right panel shows WDM 
halos. Dots are halos from the large volume run and crosses are 
from the small volume run. The dotted line shows an exact corre- 
spondence and the solid black line is a polynomial fit to the CDM 
halos. 

epochs and work equally well for both CDM and WDM 
halos. Although not shown, the errors for this fit are gen- 
erally 5 — 10%, depending on the number of halos we 
average over. Figure [TUl compares ac,wo2, from equation 
([5]), and flc from our modified form, equation ([5]), with 
dots representing halos from the large volume simula- 
tions, and crosses are from the small volume simulations. 
In the CDM cosmology (left panel), ac.yvo2 = 2. Taj-* 
with a 27% scatter. The few WDM halos with ac.wm < 1 
are also well described by this relation. In both cases, 
there are several halos with ac,wo2 > 1- These are typ- 
ically halos that first appear around a > 0.7 and have 
rapid growth phases, powered either by major mergers 
or the collapse of a sub-truncation scale perturbation in 
the WDM cosmology. These formation epochs, based on 
equation ([8]), appear unphysical since they give formation 
epochs in the period of exponential expansion when halo 
growth has stopped. Equation ^ resolves this issue by 
increasing the halo rate index, resulting in a substantially 
lowered Oc which pulls these halos significantly above the 
measured linear relation between Oc and ac.wo2- 

Figure [TT] plots the variation of Qc (top panels) and 7 
(bottom panels) with mass for CDM (left panels) and 
WDM (right panels) halos. Again, dots represent halos 
from the large volume realizations, and crosses are ha- 
los from the smaller volumes. Common to both the left 
and right panels are average trend lines for the CDM 
(solid curves) and WDM (dashed curves) halos. The 
plot shows all halos from our simulations that are well 
resolved at the end of the simulation (M200 > 400 par- 
ticles for the large volume realizations and M200 > 1000 
particles for the smaller volume realizations — see ap- 
pendix for a further discussion). At the high mass end, 
M200 > 2Mc = 2 X IO^^/i-^Mq, the average Uc's differ by 
only ~ 10%, less than the scatter for either cosmology. 
The rate index, 7, however, is about 50% higher in the 
WDM cosmologies at these high masses, reflecting the 
steeper MAHs presented in right side of Figure [H 

As mass decreases, both Qc and 7 behave differently in 
the two cosmologies. The formation epoch decreases with 
mass in the CDM runs (in accordance with "bottom up" 
structure formation) but actually increases with mass in 
the WDM cosmologies. The number of WDM halos that 
exist with A/200 Mc — all with later formation times 



TABLE 3 

RMS RESIDUALS OF MAH FITS: 



Af20o[^ ^^0] Equation ||9ll: Equation (|9)l: Equation l[8): 

g = 0.3 - 100 g = 0.3 - 1 a = 0.3 - 1 

(0.6 - 1.2) X 10^" oTi 0.044 021 
(3-6)xl0i3 0.14 0.036 0.31 
(1.5 - 3) X 10^3 ai5 0.088 0.43 



than halos with mass greater than — again suggests 
that most small halos form through an instability of re- 
gions inside larger structures. This claim is consistent 
with our picture of sub-truncation scale halos forming 
through incomplete collapse of larger perturbations. For 
CDM halos, we again see the presence of the hierarchi- 
cal structure formation from the fact that 7 is relatively 
constant throughout the entire mass range, 7 « 2 with 
modest scatter. In contrast, the WDM cosmology has 
7 increasing with lower masses, roughly as 7 oc M^q'^. 
The spray of particles with Qc and 7 much larger than 
the averages in the CDM cosmologies is a result of the 
difficulty in measuring these parameters for poorly re- 
solved halos. The MAH cannot be measured accurately 
for halos that do not grow substantially above our res- 
olution limit, and consequently our fit parameters have 
large uncertainties. The mean relation is actually within 
these uncertainties for all the low mass halos. 

We also compare our modified MAH fitting formula 
with the model of W02 in the range of halo formation 
up to the present epoch in Figure [121 In this figure, 
we compare the residuals of fits of the average MAH to 
equation ([9]), fit from a = 0.3 — 100 (dark line) and a — 
0.3 — 1 (medium line) with equation ([H]) fit from a = 
0.3 — 1 (light line) for different mass ranges taken from 
our CDM simulations. The RMS values for the residuals 
in In M are shown in Table [3l While the rms residual 
between our modified fitting formula from a = 0.3 — 100 
is approximately a factor of 2 lower than the W02 model 
for all mass ranges, equation ^ fit from a = 0.3 — 1 offers 
a substantial improvement over equation decreasing 
the rms residual by a factor of 5 or more at the plotted 
masses. The introduction of 7 is apparently an important 
correction for the early time MAH growth as well as the 
late time asymptote. Table |4] lists the values for Qc and 
7 for the profiles of Figure fT2l The fits are robust in the 
sense that Oc does not change substantially depending 
on how the MAH is calculated. In all cases, however, 
the best fit returns 7 substantially larger than 1, and 
may have some mass dependence. The exact value of 7 
depends strongly on the fit range, with fits out to a 
100 requiring a larger value in order for the MAH to 
asymptote properly. These larger values of 7 in turn 
push Qc slightly earlier. Fitting in the range 0.3 < a < 1, 
results in 7 « 1.6 with values for Oc that differ from those 
of the W02 model by only 4%. While our average halo 
MAHs all appear to have 7 > 1, there is a much larger 
spread when considering fits for individual halos. Here, 
7 ranges from 0.2 — 10 and the fits have residuals that 
are typically 20% lower than those of equation ([8]). 

4.2. Halo Substructure 



The Power Spectrum and Structure 



11 




10^2 1013 10" 10151012 1013 lOl'' 1015 

M2oo(a=100) [h-iMd 

Fig. 11. — Top Panels: The dependence of halo's formation epoch (ac in equation [9]) on mass. The left panel shows halos from the 
CDM run, and the right halos from the WDM cosmology. Dots represent halos from the large box run, and crosses are from the small 
box run. The solid and dashed curves show the trend of mass for the CDM and WDM cosmologies, and the vertical dotted line is M^. 
Bottom Panel: Same as the top panel, but now plotting the behavior of 7, the amplification factor of equation as a function of mass. 



TABLE 4 
MAH Parameters: 



M2oolh-'MQ] 


0.C.W02 


Oc 


ac 


7 


7 




0.3 - 100 


0.3 - 1 


0.3 ~ 1 


0.3 - 100 


0.3 - 1 


(0.6 - 1.2) X 10^"* 


0.54 


0.47 


0.52 


2.0 


1.5 


(3 - 6) X 10^^ 


0.53 


0.50 


0.52 


2.1 


1.6 


(1.5 - 3) X IO13 


0.48 


0.48 


0.50 


2.3 


1.8 



In this section, we compare the subhalo distribution 
of our CDM and WDM cosmologies, considering only 
halos from our smaller volume simulations. While these 
simulations do not contain a statistically large number 
of halos, the larger volume simulations do not have the 
necessary resolution to accurately describe the subhalo 
population. 

Not surprisingly, the most dramatic difference between 
the subhalo populations of our CDM and WDM halos 
is their abundance. In the CDM simulation, the aver- 
age number of subhalos with M > 1.24 x 10i°/i"iMo 
(20 particles) is roughly proportional to the mass of the 
host halo, fisubs oc M20Q at a = 1. For host halos of all 
masses, approximately 10% of the host mass is in bound 
substructures at this epoch, a value consistent with pre- 



vious studies (jKlvpin et al. |[T999f) . In contrast, for the 
WDM cosmology fisubs oc M20Q, with only about 5% of 
the host mass in bound subhalos. By a = 100 in the 
CDM cosmology, the slope of the number of subhalos 
with mass has not changed substantially, fisubs oc M2QQ, 
but many subhalos have been destroyed, and only 0.3% 
of the mass of an average halo is contained in substruc- 
ture. The steepening of the slope is caused by a more effi- 
cient destruction of subhalos in hosts with lower masses. 
Small halos today contain smaller, more weakly bound, 
subhalos than larger hosts. Additionally, these small ha- 
los undergo fewer future mergers to replenish their sub- 
halo population. By contrast, with the exception of the 
single largest WDM halo, none of the WDM halos con- 
tain any substructures at a = 100. Even this largest 
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Fig. 12. — Residuals of average MAHs to various fits. The dark 
lines are tfie residuals to equation ^ fit in the range a = 0.3 — 100. 
The medium lines also use equation but the fit is ca lculated 
over a = 0.3 — 1. The light lines fit the model of Wechsl er et al. I 
1I2OOI) . equation to the range a = 0.3 — 1. The halos are in the 
mass ranges (top to bottom) M200 = (0.6 - 1.2) X 10^"', (3 - 6) X 
10", (1.3 - 3) X lO^h-^Mg. 



halo has only fisubs = 2, as opposed to fisubs — 137 
in its CDM counterpart. Surprisingly, the shape of the 
subhalo mass function, dnsub{Msub)/dln{Msub/MHost), 
does not change substantially between today and a = 
100. The mass function is truncated at the high mass end 
sub /M Host > 0.01) and is slightly steepened. This is 
due to the increased effect of dynamical friction on more 
massive objects, pulling them towards the center of the 
halo where they are more easily disrupted and stripped 
of mass. 

Figure [T2] shows the evolution of the average number 
of subhalos in halos of various masses. CDM and WDM 
halos of all sizes show similar evolution in the average 
number of subhalos. For a < 3, mergers create substruc- 
ture, resulting in an increasing fisubsio-)- After mergers 
end, a « 3, no new subhalos are accreted and exist- 
ing subhalos gradually fall inward and are disrupted due 
to dynamical friction and tidal forces. The number of 
subhalos in the WDM cosmology, however, drops much 
more rapidly than in the CDM cosmology. The primary 
reason for this difference is the lower binding energy of 
the subhalos caused by the later formation epoch of low- 
mass WDM halos. When these halos accrete onto more 
massive halos, they are much more prone to disruption 
and consequently have a shorter life. While WDM sub- 
halos are formed with a lower binding energy/density, 
most other properties, including their average velocity, 
velocity dispersion, and mass, differ very little between 
the CDM and WDM models. The only other systematic 
difference is the average distance a subhalo lives from 
the center of its host — WDM subhalos tend to live fur- 
ther out. This is also related to the lower binding energy 
of WDM subhalos, since they are more easily disrupted 
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Fig. 13. — The evolution of the average number of subha- 
los in host halos of various masses. The solid curves indicate 
CDM halos, and the dashed WDM halos. At a = 100, the ha- 
los have masses M200 = (2 - 4) X 10^^ (light hue), (0.5 - 1.3) X 
10l''(medium hue), and > 4 X IO^^'/j.-^Mq (dark line). The dotted 
lines show fits to equation IIIOI I. 

when they move closer to center of their host. As a mea- 
sure of subhalo destruction, we have fit the evolution of 
fisub during the late-time deSitter expansion (when the 
host halos are no longer being disrupted by mergers) to 
an exponential decay, 



'^sub 



[a) 



(10) 



where the subscript m denotes the epoch where the sub- 
halo population is at its maximum, and a' is the subhalo 
decay rate, with a' = Ha = a0.0856Gyr~^h for our A- 
dominated cosmology {H — Hoo)- For our CDM halos, 
we get a = 0.38 ± 0.03, while for the WDM halos there 
is a much more rapid decay, a = 1.1 ±0.2, yielding sub- 
halo half-lives of 21 ± 2 and 7.4 ± lh~^Gyr, respectively. 
Figure [T51 makes it appear that equation (fTU)) fits better 
for lower mass halos. This, however, is an artifact of the 
fact that we have only a few high mass halos in each 
simulation, giving us poor statistics. 

4.3. Density Profile 

The radial density profile is one of the most fun- 
damental hal o properties. Previous WDM studies 
("Moore et al. ' '199 9t [Colm. Avila-Reese. k Valenzuelal 
2000; Bode, Ostri ker. fc Turok I 120011) indicate that 
WDM halo density profiles do not differ substantially 
from their CDM counterparts, and our simulations sup- 
port this finding. Figure [T?] shows average density pro- 
files for all particles bound to halos from the CDM and 
WDM cosmologies at a = 100 from our smaller volume 
runs. The solid lines are CDM halos from the mass 
ranges of Figure [51 and the dashed lines are WDM ha- 
los. The different mass ranges have been offset and the 
density multiplied by to make the differences between 
the various profiles more visible. The halos can be de- 
scribed by an NFW profile, equation(l5]) for the range 



0.05r2oo 



r2oo- Here, and ps are the the NFW 



scale radius and density. The middle panel of Figure [14] 
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Fig. 14. — Top Panel: Average density profiles times for 
all bound material in halos at a = 100 from our smaller vol- 
ume runs. Halos selected according to the mass ranges M200 = 
(2 - 4) X 10^3 (light line), (0.6 - 1.3) X lO^" (medium line), and > 
4 X IO^^/i^'^Mq (dark line), with the solid and dashed lines repre- 
senting CDM and WDM halos. The profiles have been offset from 
each other to make them easier to see. Middle Panel: Residuals 
of NFW fits to the above six profiles. Bottom Panel: Residuals of 
truncated Hernquist fits (equation 1111 ) to the above six profiles. 

shows the residuals to the NFW fits for the nine plotted 
profiles, which generall y fall in the ~ 5% range. 

Previous work (Bush a~et al. |[2005l ) has shown that well 
before a = 100, halos in a A-dominated universe develop 
clear edges at r « 4.6r2oo- Beyond r2oo, the NFW profile 
is much too shallow and the density profile for all bound 
material is better fit by a truncated Hernquist profile, 

where po and Tc are the central density and core radius, 
and rhaio is the radial extent of all bound halo material. 
Residuals to this fit are shown in the bottom panel of 
Figure [Ml This fit has larger residuals than NFW for 
< ^200: but has much smaller residuals at larger radii 
and provides a good description of the halo out to its 
actual edge. 

As with the MAHs, there is little difference between 
the density profiles of the WDM and CDM cosmologies 
for halos with M200 ^ Mc- Specifically, for the plot- 
ted profiles in this range, most differences are at the 
~ 1% level. Significant differences appear as we near 
the truncation scale, Mc, but, remarkably, the profiles 
are still well fit by both the NFW and Hernquist pro- 
files, as indicated by the residuals in the bottom panels 
of Figure [TH We define the concentration in the usual 
way, C200 = ''200 /''s- As noted above, the parameter 
rg is measured by fitting an NFW profile to our pro- 
files using logarithmically space d radial bins i n the range 
0.05r2oo — ''200 (similar to , Bullock et al. |[200lD . For equi- 
librated halos at a = 100, the steepness of the halo den- 
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Fig. 15. — Co mparison of the core radius, rc, of the Hernquist 
fit (equation 1111 ) plotted against the scale radius, Ts, of the NFW 
fit (equation [5]) from our CDM and WDM simulations. Dots 
represent halos from the large box run, and crosses are from the 
small box run. The black and dashed lines show the linear fits for 
the CDM and WDM cosmologies. 

sity profile causes to change substantially if we fit to 
a different radius. By increasing the outer fit radius to 
Tmr ~ 2.5r2oo as in Bullock et al. (2001), concentrations 
typically decrease by a factor of 2 or greater, depending 
on the mass of the halo. A more robust concentration 
could alternatively be measured as Chaio = """haio/fc from 
equation (fTT|) . although the definition using the NFW 
fit is more standard (and our new definition is not well 
defined for a ^ 3) . Figure [15] plots the relation between 
the parameters r^. and r^. The relation is well fit by a 
power law, oc with v = 0.974 ± 0.004, 0.88 ± 0.02 
and proportionality constants 0.40 and 0.24 for the CDM 
and WDM cosmologies. For what follows, however, we 
will continue to use the concentration defined with an 
NFW fit. 

Figure [16] shows the dependence of the concentration 
on mass for our halos. As with Figure [TTl the left panel 
shows CDM halos and the right panel WDM halos, with 
dots representing halos from the large volume realiza- 
tions and crosses halos from the small volume realiza- 
tions. This plot follows many of the overall trends of Fig- 
urc[TlJ Again, for A/200 ^ Mc the average values for C200 
differ by less than 10% for CDM and WDM halos and 
systematically diverge for lower masses. The low mass 
WDM halos are much "puffier" than similar mass halos 
(or their counterparts) in the CDM simulations. W02 
explained the trend of changing C200 with M200 by iden- 
tifying a relationship between C200 and Oc, which we have 
reproduced with our data in Figure [T7] Both the CDM 
(left panel) and WDM (right panel) halos follow a power 
law relation C200 oc a~^, although there is a substantial 
amount of scatter present in the CDM relation. For CDM 
halos, 13 = 0.79 ± 0.2 (sohd fine), while the WDM halos 
have P — 0.86 ± 0.3 (dashed fine), values that are both 
within 2(T from the combined slope /3 = 0.79±0.1 (dotted 
line). This finding is similar to the c (x relation pro- 
posed by W02. Note that fo r W02, r2oo is replaced wit h 
r„„ « 2.5r2oo at a = 100 l|Eke. Cole, fc Frenk I [l996l ). 
The explanation for this trend is that C200 is a reflec- 
tion of the average density at the time of collapse, so 
that halos forming earlier should have higher concentra- 
tions, exactly as observed. This trend holds even for 
WDM halos with 1/200 < Mc- This picture is further 
enforced by a toy model for concentrations proposed by 
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lEke. Navarro, fc Steinmetzl (|2001h . modified to use our 
density threshold. Noting that in WDM models C200 de- 
creases with mass below the truncation scale, Eke et al. 
(2001) defined an effective perturbation spectrum, 



d(7{M) 



(12) 



dln(M)' 

so that the spectrum also decreases at low masses. From 
this equation, a collapse epoch, ac,toy can be identified 
as the epoch where D{ac,toy)o'ef f{Ms) = l/Ca, where 
Ms = M{< 2.17rs) (the mass contained within the ra- 
dius where the circular velocity of a NEW profile reaches 
its maximum) and C„ i s a fitting parameter . A central 
density is defined as in iBullock et aTl (|2001h such that 
M200 — 4/37rrfps- We then assume that 

Ps = 200pcrit(ao)cioo (13) 
and set this density scale equal to our overdensity at 
the epoch of formation, 200pcritiac,toy)- This procedure 
yields the relation 

Ml/3 

(14) 



C200 



Pcrit (^c.toy ) 



Pcrit(ao) 

The only free parameter in this model is the constant 
C(j, which we set equal to 32. The results of this model 
are plotted as the lighter curves in Figure [TT] The model 
characterizes the simulated halos for both the CDM and 
WDM cosmologies at a = 100, and produces equally 
good agreement at earlier epochs. Our toy model dif- 
fers from that of Eke et al. (2001) only in that we de- 
fine our halos in terms of 200pcrit, as opposed to the 
epoch-dependent quantity A{a)pcrit- We find that using 
A{a)pcrit requires different values for C(^a-v) at a = 1 and 
100. 

In all of our simulations, the density profile is well 
characterized by an NEW profile for much of the halo's 
radial extent, and that the quality of the fit does not 
depend on whether a halo was taken from a CDM or 
WDM cosmology. The primary difference between halos 
in these two models — the change in concentration pa- 
rameter — follows a simple relationship with the forma- 
tion epoch that is relatively insensitive to how the halos 
form: CDM halos that have accreted their mass as viri- 
alized clumps; WDM halos with M200 > Mc that have 
formed through a period of rapid smooth accretion fol- 
lowed by growth through accretion of clumps; and WDM 
halos with Af2oo < Mc that have formed through a single 
period of rapid smooth accretion. 

5. CONCLUSIONS 

This work examines the effects of small scale struc- 
ture and merger activity on the formation and ultimate 
structure of halos in a A-dominated universe. Using N- 
body simulations with an initial power spectrum that is 
truncated on small scales, we model a WDM-like uni- 
verse where the formation of early low-mass objects is 
suppressed and compare the results to standard ACDM 
simulations. Without these seeds for hierarchical growth, 
halos above the truncation scale form through an ini- 
tial rapid accretion phase, resulting in objects with mass 
A/200 Afc = 8 X lO^^/i~^M0. These objects can then 
grow through mergers with other large halos. Many ha- 
los also form below the truncation scale through incom- 
plete collapse of larger perturbations. These halos form 



solely through a monolithic-like collapse process. Re- 
gardless of their size, WDM halos typically form later 
and faster, with a larger A'/2oo(a) than their counterparts 
in an ACDM cosmology. 

To describe this ra pid accretion, we generalize the 
IWechsler et al. I (|2002f ) mass accretion history formula 
by introducing a rate index, 7, that controls the growth 
rate evolution, dln(M)/(iln(a) oc a""*". This parameter 
is necessary in order to fit the halo MAH into the far fu- 
ture. Otherwise, the halo mass approaches its asymptote 
much too slowly. Additionally, for ensemble-averaged 
halo histories fit to the present epoch (a < 1), we recover 
7 ~ 1.6, substantially higher than the 7 = 1 assumed in 
IWechsler et aTl (|2002D . The lar ger 7 improves the mean 
MAH fit by a factor of 5. Our generalization also reduces 
the number of objects that have unphysical formations 
epochs, Oc ^ 1. 

We have also calculated abundances and decay rates 
for substructure in our dark matter halos. As expected, 
WDM halos contain much less substructure, and most 
host halos destroy all of their subhalos by the end of our 
simulations at a = 100. The decay of subhalos is (on 
average) exponential in time with half-lives 21Gyr and 
7.4Gyr for the CDM and WDM cosmologies. This is 
consistent with the picture discussed below where later 
formation times create lower concentration objects that 
are more prone to disruption. 

Despite differences in the formation process and sub- 
structure abundance, WDM halos exhibit NEW den- 
sity profiles, just like CDM halos, albeit with a differ- 
ent C2oo(-^^20o) relationship. Halos in WDM cosmologies 
have lower concentrations than their CDM counterparts, 
but follow the same 0200(0,0) relation, allowing us to relate 
concentration to mass in both cosmologies using a 1 pa- 
rameter toy model that characterizes the concentration 
using the linear power spectrum. This characterization is 
motivated by the idea that the concentration is set by the 
cosmological background density at the epoch of collapse. 
The form of the halo density profile persists, even though 
the method of formation and amount of substructure is 
changed substantially in our WDM cosmology. Taken 
together, our results suggest that the form of the halo 
density profile is set not by merger activity, but instead 
through large-scale modes of the gravitational relaxation 
process. Halos appear to be very efficient at erasing their 
initial conditions and do not care whether their mass was 
accreted rapidly, slowly, in clumps, or continuously. The 
only aspect of the density profile directly linked to the 
halo formation process appears to be the concentration, 
which set by the formation epoch, the epoch when the 
mass accretion rate drops to a specific value. In particu- 
lar, substructures and major mergers seem to have little 
effect in driving the equilibrium structure of a halo. 
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Fig. 16. — The concentration (c200 = fs/»"20o) as a function of mass for CDM and WDM halos. Dots represent halos from the large box 
run, and crosses are from the small box run. The dark solid and dashed lines show average C2oo(A^20o) values for CDM and WDM halos, 
and the vertical dotted line marks the truncation scale, Afc. The light solid and dashed lines show the calculated C20o(JW20o) relation from 
the toy model discussed in §4.3. 
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Fig. 17. — The concentration as a function of formation epoch, from equation for CDM and WDM halos. Dots represent halos 
from the large box run, and crosses are from the small box run. The solid line is the trend-line for CDM halos, and the dashed line is the 
trend-line for WDM halos. The dotted line is the trend-line for the combined sample of all CDM and WDM halos. 
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Fig. 18. — Left Panel: The fraction of mass in a collapsed structure of a given mass as a function of mass for a number of WDM simulations 
with varying degrees of mass particle resolution (A/pfj^f^^^g = 6.20 X 10** - 3.97 X 1010/i-1Mq) and simulated volume (50 - 200/i-iMpc). 
The vertical dotted line represents the truncation scale, Mc- Right Panel: The same mass fractions, now plotted mass in units of number 
of particles. 
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APPENDIX 
MASS FUNCTION 

The upturn in the FOF mass function (Figure [6]) at small mass is a surprising and potentially troubling feature of 
our WDM simulations, and it is necessary to determine if this phenomenon is a result of the simulation or of the FOF 
group-finding algorithm. We have conducted a series of 5 WDM simulations with various volumes and mass resolutions 
to examine this effect, and present the resultant mass fractions at a = 1 in Figure [181 Here, we are plotting the mass 
fractions, as opposed to the mass functions shown in Figured The mass fraction and mass function are related by 
dn{M) / d\n{M) = {p/AI)f{AI). We used a mass resolution of 3.97 x 10^"h~^MQ as our base model, conducting three 
runs at this scale with comoving box lengths of 200, 100, and 50 /i^^Mpc (solid, dash-dot, and dash-dot-dot-dot lines), 
and two additional runs with box length 50/i~^Mpc with mass resolutions 6.20 x 10^ and 4.96 x 10^h~^MQ (dotted 
and dashed lines). The left panel shows the mass fraction as a function of mass in units of /i~^Mq, while the right 
panel shows the same fractions in units of number of particles. 

The most immediate feature of the left panel is that the location of the upturn changes with mass resolution, 
independent of the simulated volume. The right panel of Figure [18] further indicates that this upturn is a purely 
numerical artifact. Here, the mass fraction is plotted in units of number of particles in the halo, and we see that at 
low masses all resolutions converge on the upturn. Finally, we also note that for the two runs that have corresponding 
CDM cosmologies (solid and dotted lines) the location of the minimum in the mass fraction [Mmin) corresponds to a 
transition point: halos with M > M^^^ have corresponding halos in the CDM simulation, while those with AI < Al^m 
do not (see Figure[7|). In these, 90% of all halos with A/IpoF ^ Mmin have corresponding halos in the CDM cosmology, 
while less than 1% of halos with AdpoF < Mmin have identified cross halos. Even though the input power spectrum 
for the CDM and WDM cosmologies is vastly different, it is not the value of Mc, but rather the of value Admin that 
determines the mass where we are unable to find corresponding halos in the CDM cosmology. This trend is especially 
noteworthy in the case of the small volume realization, where Mmin is almost two orders of magnitude below Ale- 
Finally, the upturn appears substantially less pronounced for the SO mass function. Most of the FOF halos below the 
upturn have no central overdensity greater than 2QQpcrit but instead exist as a more linear string of particles that are 
not even identified as energetically bound by subfind . 

One expects an upturn eventually, since a FOF group finder puts every particle in a group of some size. At any 
given time, 50-70% of all particles end up in "groups" of one particle, requiring an upturn somewhere, something that 
you would expect even in a standard ACDM simulation (although the effect there should be much smaller due to the 
large number of actual groups at this mass range). Our results suggest that this contamination is an issue for groups 
containing fewer than 

(lifer) "°" (Ai) 
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particles when a hnking length of 0.164 is used. 
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